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Abstract —We propose optimized 2D Non-Uniform FFT 
(NUFFT) algorithms and apply them to Magnetic Resonance 
Imaging under radial acquisitions. 

Index Terms —NUFFT, MRI, Filtered Backprojection 

I. Introduction 

In MRI, much effort has been devoted to acquisition meth¬ 
ods using non-Cartesian grids [1], such as radial [2] and spiral 
[3] scannings, to achieve fast data acquisition. Spiral scans are 
motivated by the simpler generation of field gradients which 
do not require high-speed, abrupt changes [3] as compared 
to the Cartesian one. Radial scans are also less prone to 
residual motion and changes of the imaging plane during 
interactive scanning making them attractive tools for image- 
guided interventions [4]. 

However, the use of standard FFTs to process non-Cartesian 
data breaks down under radial or spiral scans so that 
MRI imaging under alternative collection procedures requires 
more involved data processing to properly balance accuracy 
and computational burden. Fortunately, Non-Uniform FFT 
(NUFFT) provide a convenient trade-off between accuracy and 
complexity [1] [5] [6] [7] [8]. 

Recently, an optimized approach, based on a general and 
new perspective, has been presented in [8] for all the three 
types of NUFFTs. The choice of the window function has 
been optimized to obtain more accurate results than those 
available in the literature. The computational costs and the 
memory requirements of the proposed schemes have been 
theoretically analyzed and their very convenient performance 
assessed, theoretically and numerically, by comparisons with 
a number of relevant and popular approaches. The proposed 
method has proved to be more accurate than the confronted 
schemes without burdening the computational and memory 
requirements. However, the approach has been limited to ID. 
The purpose of this paper is twofold: 

• extend the approach presented in [8] to 2D; 

• apply the “optimized” 2D NUFFT to non-Cartesian k- 
space acquisitions; in particular, radial MRI processing 
is considered. 

II. NUFFT-based MRI Reconstruction 

The signal r acquired by the coils of an MRI system is 
related to the spin density p by the following relation 


r(k x (t),ky(t)) = JJ S(x,y)p(x,y)e~^ k ^ x+k y^dO, 

° ( 1 ) 
where S is the receiving coil sensitivity map, k(t) = 
(k x (t),k y (t)) denotes the k- space trajectory and O is the 
Region of Interest (ROI), having considered, for the sake of 
simplicity, a 2D case. 

Using a pixel-based approximation for p and sampling the 
signal r at the time instants t n , eq. (1) is regarded as a system 
of linear equations 

Ap = r, (2) 

where r is the data vector whose generic element is r n = 
r(kx(tn),ky(t n )) 9 p is the unknowns vector containing the 
pixel values p s = p(x s ,y s ), with (x s ,y s ) the pixel centers. 
We assume that the pixel centers (x s ,y s ) are arranged on a 
Cartesian grid, namely, (x s ,y s ) = (x p ,y q ), with s = s(p,q). 
The generic element A ns of A is provided by 

A ns = S(x s , y s )e-^ kx ^ x °+ k y( tn ^AxAy, (3) 

with Ax Ay the pixel size. 

A pseudosolution to the system (2) is typically searched for 
by solving the associated Euler equation 

A} A p = b (4) 

where b = A} r and A ^ is the adjoint matrix of A whose 

generic element is 

4 „ = S*(x s , AxAy, (5) 

where * denotes complex conjugation. Eq. (4) is then solved 
by using iterative techniques minimizing the following error 
functional 

*(£) = \\A'Ap-A'tW 2 , ( 6 ) 

for example, the conjugate gradient method. 

The matrix-vector multiplications involving matrices A and A ^ 
can be computed by 2D NUFFTs of Non-Equispaced Results 



(NER) and Non-Equispaced Data (NED), respectively. We 
now briefly recall the NER and NED Non-Uniform DFTs in 
the ID case, for the case of simplicity. However, the results 
in the next Section will regard the 2D case. 

In the ID case, the NER NUDFT of uniform samples 
{ z k}k= 2 -^N/ 2 ’ eva l ua ted a t non-equispaced grid points xi E 
y], is defined as 

N_1 

zi= z k e~* 2 ™'*, / = M. (7) 

fc=-f 

On the other side, the ID NED NUDFT of samples {zi}f£ l9 
located at non-equispaced points and evaluated on an 

equispaced grid, is defined as 


M 

z k = Y^ zie~ j2 ™ 1 &,k= —N/2, ...,N/2-l. (8) 

1 = 1 

III. Numerical Results 


table II 

RMS AND MAXIMUM ERRORS: DOUBLE PRECISION ARITHMETICS. 



Type-2 NUFFT Algorithm 

KB 

Optimized 

RMS 

5.40 • 10“ 1U 

8.87- lO -11 

MAX 

7.92 • 10- iU 

1.73 • 10- iU 



k x t" 1 ' 1 ] 


In the present Section, we focus the attention on the exten¬ 
sions of the NER- and NED-NUFFT algorithms developed in 
[8] to the 2D case. 

The accuracy of the NER-NUFFT is quantitatively assessed in 
the case of 12 x 12 non-uniformly and randomly distributed 
sampling points generated in the [—6,6]m _1 x [—6,6]m _1 
spectral region. The data to be transformed have been assumed 
complex, with real and imaginary parts uniformly distributed 
in [—0.5, 0.5]. Tabs. I and II illustrate the accuracy in terms 
of Root Mean Square (RMS) errors and maximum errors 
as compared to the respective NER-NUDFT for single and 
double precision arithmetics, respectively. In both the cases, 
the accuracy of the optimized NUFFT is compared to the 
scheme in [5]. For single precision, the optimized NUFFT 
is compared also with that in [1] which is available for single 
precision only. As it can be seen, our approach outperforms 
the compared ones. 

The performance of the NED-NUFFT is illustrated by showing 
the result of the computation of Ap which amounts at the 
so-called Filtered Backprojection (FPB) algorithm for radial 
MRI acquisitions. Fig. 1 illustrates the radial spectral sample 
locations, while Fig. 2 depicts the reconstruction of the Shepp- 
Fogan phantom. Accuracies similar to those in Table I are 
achieved. 

TABLE I 

RMS AND MAXIMUM ERRORS: SINGLE PRECISION ARITHMETICS. 



Type-2 NUFFT Algorithm 

KB 

Optimized 

Fessler 

RMS 

6.23 • 10 -4 

1.11 ■ 10 -4 

122■10“ 4 

MAX 

6.93 • lO -5 

1.12 ■ 10“ b 

140- 10“ b 
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Fig. 2. Phantom reconstruction. 
















